(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     51117,       1218]
NotebookOptionsPosition[     49706,       1186]
NotebookOutlinePosition[     50062,       1202]
CellTagsIndexPosition[     50019,       1199]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{"ClearAll", "[", "\"\<Global`*\>\"", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Theta]", "=", 
    RowBox[{"\[Pi]", "*", 
     RowBox[{"(", 
      RowBox[{"7", "/", "180"}], ")"}]}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "Canting", " ", "angle", " ", "of", " ", "the", " ", "applied", " ", 
    "field", " ", "with", " ", "respect", " ", "to", " ", "z"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Epsilon]", "=", "0.028"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "Microwave", " ", "power", " ", "is", " ", "taken", " ", "to", " ", "be", 
    " ", "much", " ", "larger", " ", "than", " ", "its", " ", "actual", " ", 
    "value", " ", "just", " ", "for", " ", "illustrative", " ", "purposes"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]", "=", 
    RowBox[{"2", "\[Pi]", "*", "0.395"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"Fixed", " ", "driving", " ", "frequency", " ", "395", "GHz"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Phi]", "=", "0"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Right", "-", 
     RowBox[{"handed", " ", "polarization"}]}], ";", " ", 
    RowBox[{"\[Phi]", "=", 
     RowBox[{
      RowBox[{"\[Pi]", "/", "2"}], " ", "for", " ", "linear", " ", 
      "polarization"}]}], ";", " ", 
    RowBox[{"\[Phi]", "=", 
     RowBox[{
      RowBox[{"\[Pi]", " ", "for", " ", "left"}], "-", 
      RowBox[{"handed", " ", "polarization"}]}]}]}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Delta]m", "=", "0.001"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "Initial", " ", "spin", " ", "fluctuation", " ", "to", " ", "reduce", " ", 
    "incubation", " ", "time"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Omega]E", "=", 
   RowBox[{"2", "\[Pi]", "*", "1.317"}]}], ";", 
  RowBox[{"\[Omega]K", "=", 
   RowBox[{
    RowBox[{"-", "2"}], "\[Pi]", "*", "0.026"}]}], ";", 
  RowBox[{"\[Omega]a", "=", 
   RowBox[{
    RowBox[{"-", "2"}], "\[Pi]", "*", "0.001"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Alpha]", "=", "0.005"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Gilbert", " ", "damping", " ", "is", " ", "amplified", " ", "to", " ", 
     "achieve", " ", "the", " ", "steady"}], "-", 
    RowBox[{
    "state", " ", "oscillation", " ", "quickly", " ", "and", " ", "reduce", 
     " ", "computation", " ", "time"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"tmax", "=", "2000"}], ";", " ", 
  RowBox[{"(*", 
   RowBox[{"Duration", " ", "of", " ", "simulation", " ", "2", "ns"}], "*)"}],
   "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"h", "[", "t_", "]"}], ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"h0", " ", 
       RowBox[{"Sin", "[", "\[Theta]", "]"}]}], "+", 
      RowBox[{"\[Epsilon]", " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Omega]", " ", "t"}], "]"}]}]}], ",", 
     RowBox[{"\[Epsilon]", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"\[Omega]", " ", "t"}], "+", "\[Phi]"}], "]"}]}], ",", 
     RowBox[{"h0", " ", 
      RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"mA", "[", "t_", "]"}], ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"mAx", "[", "t", "]"}], ",", 
     RowBox[{"mAy", "[", "t", "]"}], ",", 
     RowBox[{"mAz", "[", "t", "]"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"mB", "[", "t_", "]"}], ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"mBx", "[", "t", "]"}], ",", 
     RowBox[{"mBy", "[", "t", "]"}], ",", 
     RowBox[{"mBz", "[", "t", "]"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LLG", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mA", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"h", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mA", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mB", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mAz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mAx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mA", "'"}], "[", "t", "]"}]}]}]}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mB", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"h", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mB", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mA", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mBz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mBx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mB", "'"}], "[", "t", "]"}]}]}]}]}], "]"}]}], "}"}], "//",
       "Flatten"}], ")"}], "~", "Join", "~", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"mAx", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mAy", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mAz", "[", "0", "]"}], "\[Equal]", 
       SqrtBox[
        RowBox[{"1", "-", 
         RowBox[{"\[Delta]m", "^", "2"}]}]]}], ",", 
      RowBox[{
       RowBox[{"mBx", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mBy", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mBz", "[", "0", "]"}], "\[Equal]", 
       RowBox[{"-", 
        SqrtBox[
         RowBox[{"1", "-", 
          RowBox[{"\[Delta]m", "^", "2"}]}]]}]}]}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"variables", "=", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"mA", "[", "t", "]"}], ",", 
      RowBox[{"mB", "[", "t", "]"}]}], "}"}], "//", "Flatten"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"res", "=", 
   RowBox[{"ParallelTable", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"h0", "=", 
       RowBox[{"2", "\[Pi]", "*", "0.001", "*", 
        RowBox[{"(", 
         RowBox[{"j", "-", "300"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"sol", "=", 
       RowBox[{"NDSolveValue", "[", 
        RowBox[{"LLG", ",", "variables", ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0", ",", "tmax"}], "}"}], ",", 
         RowBox[{"MaxStepSize", "\[Rule]", "0.001"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Is", "=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "2", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "1", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "1", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "2", "]"}], "]"}], ",", "t"}], "]"}]}]}], ")"}], 
        "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "5", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "4", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "4", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "5", "]"}], "]"}], ",", "t"}], "]"}]}]}], 
         ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"NIntegrate", "[", 
       RowBox[{"Is", ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", 
          RowBox[{"tmax", "/", "2"}], ",", "tmax"}], "}"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "0", ",", "600"}], "}"}]}], "]"}]}], ";"}]}], "Input",\

 CellChangeTimes->{{3.784226962268992*^9, 3.784227228761891*^9}, {
   3.7842275027982073`*^9, 3.784227589442974*^9}, {3.784227655152667*^9, 
   3.7842276554218817`*^9}, 3.784227851779511*^9, {3.784228109319489*^9, 
   3.784228151830392*^9}, {3.784228353040759*^9, 3.784228360587185*^9}, {
   3.784228438570054*^9, 3.78422843886664*^9}, {3.784228529783247*^9, 
   3.784228556552698*^9}, {3.784228614972218*^9, 3.784228707442814*^9}, {
   3.784228772223241*^9, 3.7842288038547077`*^9}, 3.7842288833757763`*^9, {
   3.784228964278302*^9, 3.7842290547253933`*^9}, {3.784229105687675*^9, 
   3.784229211003819*^9}, {3.784229301049076*^9, 3.784229538772077*^9}, {
   3.784233261327928*^9, 3.784233290456541*^9}, {3.784233509907217*^9, 
   3.78423351924338*^9}, {3.784233643898909*^9, 3.784233644704205*^9}, {
   3.7842338436015463`*^9, 3.78423384415106*^9}, {3.7842338859399776`*^9, 
   3.784234050734558*^9}, {3.784234117964368*^9, 3.784234220052305*^9}, {
   3.7842346088995457`*^9, 3.784234623687046*^9}, {3.784234664690909*^9, 
   3.784234701419993*^9}, {3.784234898411098*^9, 3.7842349757371283`*^9}, {
   3.7842350472771473`*^9, 3.784235061835257*^9}, {3.784235166316753*^9, 
   3.7842351732323503`*^9}, {3.784235251324011*^9, 3.7842352574040956`*^9}, {
   3.784235289484105*^9, 3.784235327281433*^9}, {3.784235375854424*^9, 
   3.784235376396117*^9}, {3.784235649852745*^9, 3.784235650626583*^9}, {
   3.78423569553428*^9, 3.7842356957866993`*^9}, {3.784236015894725*^9, 
   3.7842360161375713`*^9}, {3.784236053732073*^9, 3.784236101939043*^9}, {
   3.784236156334631*^9, 3.784236207295868*^9}, {3.784236285437698*^9, 
   3.7842362866907177`*^9}, {3.784236327381797*^9, 3.784236385351768*^9}, {
   3.784236520622822*^9, 3.784236541128675*^9}, {3.784236619922337*^9, 
   3.7842366585641336`*^9}, {3.784236759451808*^9, 3.7842368020890417`*^9}, {
   3.7842368533120537`*^9, 3.784236908228676*^9}, {3.784236950202259*^9, 
   3.784237210409375*^9}, {3.784237330445344*^9, 3.7842373643425093`*^9}, {
   3.784242346085524*^9, 3.784242700084516*^9}, 3.784242803084869*^9, {
   3.7842428335804234`*^9, 3.7842429554282722`*^9}, {3.784487858736527*^9, 
   3.7844879036983147`*^9}, {3.784487938035633*^9, 3.784488013832548*^9}, {
   3.784488079260935*^9, 3.784488112647812*^9}, {3.784488232448867*^9, 
   3.784488259390422*^9}, {3.784488307939101*^9, 3.7844884442069387`*^9}, {
   3.784488591758436*^9, 3.784488591946334*^9}, {3.7844886250955133`*^9, 
   3.784488661072905*^9}, 3.784488693129335*^9, {3.784488745200171*^9, 
   3.784488745533964*^9}, 3.784489223120585*^9, {3.78451312066894*^9, 
   3.784513213646172*^9}, {3.7845132715992002`*^9, 3.784513272524568*^9}, 
   3.7845133378267393`*^9, {3.784513472676453*^9, 3.784513477048867*^9}, 
   3.7845769620912733`*^9, {3.784577088884696*^9, 3.7845770969403677`*^9}, {
   3.784577134205518*^9, 3.784577185778541*^9}, {3.7845772229635477`*^9, 
   3.784577261196495*^9}, {3.784577291750091*^9, 3.784577388440825*^9}, 
   3.784577614211357*^9, {3.7846574629201317`*^9, 3.784657463387754*^9}, 
   3.784657500864352*^9, {3.784657532424686*^9, 3.784657540759239*^9}, {
   3.784657594267646*^9, 3.784657617996209*^9}, {3.784657700354312*^9, 
   3.784657744489027*^9}, 3.784657785603743*^9, {3.7846578371437817`*^9, 
   3.7846578377700567`*^9}, {3.784657888400833*^9, 3.784657888817245*^9}, {
   3.7846579255207663`*^9, 3.784657973082341*^9}, {3.784658006523283*^9, 
   3.784658007889015*^9}, 3.784658054151285*^9, {3.784658086438966*^9, 
   3.7846580870994997`*^9}, {3.7846583908921537`*^9, 
   3.7846583918218536`*^9}, {3.7846584357558284`*^9, 
   3.7846584361137543`*^9}, {3.784658741998811*^9, 3.784658849602047*^9}, 
   3.7846588806535096`*^9, {3.784658930815446*^9, 3.7846590110481186`*^9}, {
   3.7846590452214074`*^9, 3.784659048077553*^9}, {3.7846605008772554`*^9, 
   3.7846605089242964`*^9}, 3.7846619574817944`*^9, {3.7846619962482967`*^9, 
   3.784661996888412*^9}, {3.784665048917243*^9, 3.784665059385776*^9}, {
   3.78478978246836*^9, 3.7847897924477453`*^9}, {3.784789835490391*^9, 
   3.7847898464330835`*^9}, {3.7848380445901365`*^9, 3.784838046855751*^9}, 
   3.784841176388877*^9, {3.791128131900264*^9, 3.791128299082412*^9}, {
   3.7911283540485783`*^9, 3.7911283545589314`*^9}, {3.791128504055644*^9, 
   3.791128507873851*^9}},ExpressionUUID->"5981de4b-bc46-466a-bea8-\
49eb3958bb2a"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"res", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", "Blue"}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", 
   RowBox[{"FrameStyle", "\[Rule]", "Medium"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7846767898799896`*^9, 3.7846768101199427`*^9}},
 CellLabel->"In[52]:=",ExpressionUUID->"a729b156-30d0-405c-ba5f-88b108eccbfe"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {RGBColor[0, 0, 1], PointSize[0.0055000000000000005`], AbsoluteThickness[
     1.6], LineBox[CompressedData["
1:eJw12XlcjOv7B/ApHJEllEKYFmtUlohKnyiUykhplWnfa9prpuWZQWVJZUmH
MBEJJfsWhkLKEml1MCRa0OZwEH5+367xB6+3u6trPs9yP/dzj4ZnmJ2PPIvF
cv7z1///2/+ny5TNVc3H1v036T+QtNB0nr1JLlkBRjuPVRh9kVkJm1xrUzUu
HCArI6ZpabwO/yBZDRsdrp7hWxwiq6Ox7uLP3YpiMhtyfCW7o7EyayJ/3rNf
Pe9k1sZtY/VDcc555Kloq6oJ1Xks83RkPWvR7l1xmDwTcasQerJM5lmYNmrA
yolLj5B1UV2sGilXJrM+Gt+kOb2yyCfPQaPk+xuzSpnnQqA5PyV19VHyPMT2
jm8qeibzfOivT5qR6nSMbABW3+P6R40yL8Duw4PvV64rIC/Egl07slMfyGyI
/JGqtUmGx8mLEO68Z0F+jsyLcfvawhFj2mU2QvvtWRlbxxeSjeHm72Jcriez
CUYpRF6JmizzEmzcwb74qVZWb4pJKnYcfyuZAc8M7sVTgfR5GGBog6koSlH2
+c2goNUeN2EaHQ/GDLeH8rp2hsuO/1Lojvk+sHI9XS/MUrQtrXU+GJlF48ug
PfZSYre1P40vg5H8HYM2yx2m/ePm+LqhO21jfX6/GXMcdlqm01x1jsYtYDEq
2lvfTULjFgjJ4O7TyblP48sRZHF0v7i0hsaX4+f+MZWTnV/S+ApMX1vYdYjz
nsZX4EPpPMtFN7ppfCUsC31Uf175RuMrsbrBfv9MIxb6xy2RqjrPIiLpN41b
wr9dqbam5QvVW+F+zaENVm4faNwK5fyVC54lSGl8Fdas+x0R9O0Jja+Civwe
+bPBlIdljc/KLS5o2kfj1nCc9Ldi80TZ8beBRUOumkHqZTp+NvgwaV/0ec4t
GreFrsIVY5baHRq3RZ+c26Bq7woaXw2jz29zTR9X0vhqcH9sT505+SGNcxDd
mlp7oe5Rv8FBa9ZAsznLq+nnOfioejpsYw9ZwoHmHnak+54nVL8GRqqbOXI6
T6l+Dbza9ZfFHyMza6Cip6V0fEoN1a+B+Na1oFm7ySw7rLqRomnylQw7VA4J
unza9RnV26HKYkX65YtkiR0G6uWukx9SS/Vr8eBnEueGCxlrwXRL9USFZGYt
tEcdWtbcTZashXfJibZ7C+qo3h722tODS/hk2GONV0BSWCmZsUdgZtvN6u9k
iT3elmr/VbCgnuodIBX7Rc0KJ8MB81Ier884QWYcYFqalb7+DVnigNG+Ow/y
1Bqofh1aTvF3T7AhYx1m1E3ZcYchM+ug3Nx0c/85smQdhudw7HNaZPWO2O1Q
N6hkbCPVO8LVLnP8p+VkxhHioXV/hceQJY54oZC3zPkomeWEWxsLuNdrZPVO
kFyMVy+Ua6J6JwieDmVp6ZElTtjWM+TXHDcyyxmtgQsevUojwxlv94YZG16Q
1Tsj0vf18hWvZfXOGKd6J1d1+HOqd0FWno/qWUMyXBBR1cea5E1mXFCW5n/S
I4MscYFJUMc30VVZvSueqE0oTG2R1bui6CzXLVrpH6p3heCzTtoaI7LEFctP
Zq/V8CWz3BBxokbUlkmGG2bZdQYXX5PVu6H39i9W1DtZvRu6r3Q+MBn1gurX
41b1lyEjjclYjw+zbPLbfcnMenwxTkNNFlmyHnlcP5OqUlm9Ozqq8s49ey+r
d4dl8qqI7tEvqd4dNsrpfZpLyBJ3cHJ6NgUGkFkbwJh3tFbsJmMDfNZvPblE
IqvfgEKJZ2dVh6x+A7aYqYXxVF9RPRdxidOjdZaR2VxcNSiv/hFKBhd2I/Tv
v9pH5nJh7MIfWXuXzHCxydL6VEMPWczFqWs50ZEjpdSPC2++Wl7iFLKUC6XX
+8enGZFZHphnd2pfxhoy2wPLVXONM/3I8EDAoav+mxPJXA9ocJdfDdhFZjxg
d4zzZkEhWeyB7Y+UU1puyPp7YGKJeFzIM1l/D8zVG3WprE3W3xP/PRf1vPgl
6+8J9ZeuRcVjXlN/T5zXGPFh6nQy1xPOS5afNzAmM554t8rk573VZLEnNGNy
blV4kiWeuJLGt9KIJks98U3Jrv5KCpnlhTE/dtfy95LZXmjKm99pWyDr74UV
SxNWTb0o6++Fjex5IT/KZP29cDx74+imall/L0wLi1/W+o+svxd2SY2nhbfK
+nsheR97fFevrL83cp9oVKmx3lB/b8jr7bCoH06GNxQFzBuzSWSuN+ZFjdKM
n0tmvMF//G9m7Cqy2BtBD+x1XwSQJd4oNzHZWJFOlnpjwJvi8w8vklk+cAjf
Vzfyray/Dw4Xj2sKV26m/j6Y+uz80gfLyVwfTNuc5aHKJzM++L5xSJWoiCz2
QXJ2rqtASpb44LjVX5lHR7+l/j54abbe9ekyMssXn5nkrz2RZLYvTusN1hpx
hAxfbI9bdHHhEzLXF+VZa7O2/iIzvphkZZSirdNC/X0xzJAr0XIkS3wRqHjY
+KiILPUFNOu/lhWRWX6Y2nJlZ04Dme0HT433JosGvKP+fhj9K+Pxldlkrh+i
vO07dZzIjB/ym3+vOiQii/1wd0VdjUYRWeKHbw3GHTfryVI/rGnceGiT/Hvq
748Hsyufxc8ms/1x/mJj6HEnMvxR/qglSm0TmeuPsEVl+k9Okxl/9FaqWjU+
J4v9YZ239sUChVbq74+2XS6ZnfPJUn/8W6h4RMGTzApAX8SvXykZZHYAEp62
Lgu7TkYAiv+reXqng8wNgNfWkvod49uofwBStp9c88iSLA4Ar8RtfFo8WRKA
1xzlDzcLydIAJG529ItvIrMC8Tk9X/WyYjv1D4S/eXVssgkZgah6u/3fp2Fk
biDGZtd9vnCYzARC/0rzU506sjgQ86+Xti4c2kH9A5H47Vj6iyVkaSCSHKeM
Y0eRWUEwM1VpkT9BZgfBrdZh9RYpGUE44xN87aLqB+ofhGFGrN07VpOZILx/
vMJaOY0sDsLxMCuPVbfIkiD4XNd2XfiDLP3z+xW0Hz03+Ej9g2GZ1upqHE5m
B6Ong3ferYiMYEwRaO1FO5kbjNjUmzUd0z5R/2AUWLiv8vQli4NxyXN9d8FR
siQY4VfsH0hayNJg6Dr2fimZ2kn9QzB/zTdxgj+ZHYK5d0YrTT9JRghWez4W
X/pE5obgZceN5r3juqh/CMJT9qRZGZPFIRgRajI6ZgNZEoIBm5xfq28kS0Ng
+2O1t1kBmRWKf54/UH1XRWaHIqfdNHV4NxmhKJjtZXN1bDf1D8WHjdeVO43J
TCjMG6blFHuRxaGwjwk98mMrWRKKsbEVYxrPkqWhEBYnNFg+J7PCMKg7vpcz
sIf6h0GS0B3fPZuMMEiH/dqt60TmhqFvuJ1IcSOZCUPvp8bd24rJ4jC8HsU3
ONtEloRBy+Fes2hwL/UPg2F65QI5AzKLh9lby17oe5GVeNA+9VusspPM5mGQ
9QnPc7fI+jz43iucPKKHDB6yGvd0zdD63G8OD/sq9gxVdCBzeZiR5y65kErm
8bBTxzVP9xqZ4WHi1dGP4zvJmTzovdq2O1f7X8rDg/mr4f77XcglPBg8UHUQ
ZJElPDy05saY3SdX83CEpz2sV+4L5eXhv6LYxD1G5C4etiXlOevGkFnh2Htz
r3PZGbJSOF4aRHq4fCKzw1Fx8bfwm85Xyh8OS1fP78cCyQjHfKu7iwJOkDnh
SPBVfb28g8wNx7CCzvLls/+j/OEofOLfGcYjM+GIDtU9fP88OTMcfz9Mc3L9
ThaH49vtnPuzzL5R/nBIp2W1W2whS8JRV83uLn5Krg7HkS2fHHgTv1P+cKz9
R6t+ZwC5Kxwt9uaCyZfIrAhETpHumvzXD8ofgdgTJyS568jsCJz6sLlHfJys
H4E7j++GLugjIwKGY48+9VvTR/kjcDL07lgcJ3MjYKT7Y2cl6yflj4C9TeH3
0a5kJgK4G5wy5RI5MwI/q8uLRqj8ovwRmMgWRj+PIpdEwOSaMefvOrIkAn36
+QYOi39T/gjo/trsqC4mSyNQFOtkOXAgS9KfPwKLvi27M2wRmRUJm6kjjs0L
IytFYt2qR2uFx8jsSOT09dp2vyTrR2Lxs6GXzqjJ9RuRsPX70XZsLZkTCd8v
wfXvMsjcSFRfSEoRPCTz/oxLiozDhsn3m4lEXZFj+wNrcmYkDOQ3nTq0gyyO
xNyPnrq9T8glkTBX1NpeM3ZAvyWRaOfKu9usJ1dHwvW03Y2Qo2RpJKTOTauX
dpK7ImF6O928cfFAyh+Fy13nvM3SyEpRKH0RMHZrPZkdhaUNhT5l0wdR/ij8
Pn552C8BGVEI0NdMXPWEzInCWjvjtuvT/6L8UZiw3dzWX0TmRaFda1iy+0sy
E4UYe3+9C8aDKX8UGj4tZokOkMVReDNx/fbnLAXKH4XXEctiGn3Jkiio3s2O
2f2YXB2FI+UngqYaDaH8UdD8d2TgiUJyVxRseNMeGU8YSvmjEXmp2/ZHBlkp
GkH/qI8epKBI+aPRNiHuFn8jWT8aozvv5yXJD6P8f35e0/7k/M1kTjSe9g3P
v604nPJH4+Cdv4042WReNIJ3RMwYPGUE5Y+GSUif9ZjL5MxoDGAbVh9ePZLy
R+P1mAfW7zrIJdGIedQbn7ZcifJHY9iAdLdTf5Oro6G8zxCCTrI0GvZQmT5m
5SjKHw31vTVGe46QWTGYGHgjbZL8aMofA/2PWxvqvMjsGOy/ZTrhTQVZPwbJ
X8Pe+80ZQ/ljkDfi0vSdB8icGLgFu1ilD1em/DGYcsLdJ01I5sVAPScr4sI3
MhODck4zYxmjQvlj0LTJpyvhK1kcgy2+SWe3JY6l/DFQG16VcUFBlfLHYEiX
ddSSv8nVMSiqjr0dqqdG+WNwx7S3aFcVuSsGmr0XDv4OGUf5Y+H4M/vAUNXx
lD8W4/5Wrxp9j8yOhWLx89cJSRMofyykc3bcOW+iTvn/WHE9W3fQRMofC/cR
30f61JK5sQgWrnj7sngS5Y9F14kz3TN3Tqb8sVBefau025RN+WPRURns4FhK
FsdCyT3l1aElGpQ/FrUVKeNn3SNLYlF+XbrXe50m5Y9FZUlu4tUPZGksWoeO
33B8qxblj4XWgjO65/S1KX8c1F7qWge9IivFoSxY63TI3imUPw4eGwWmPs5T
KX8ctH3zLIdNnUb54xD+vlBH9xeZE4ejtybMK3wznfLHYWXyi8u8ZzMofxw0
DTuEGc9mUv442EpPhM18q0P541CUrbFt48DZlD8OCmsOxv1coEv545Dvsshr
cqIe5Y+Ddcb8K9bB+pQ/Dqr1czvzp8yh/HG4LvyxPq+D3BWHt0oSG/+bcyl/
PJTStUeYHZ5H+eMh1/HROmfPfMofj24rn5GD9htQ/nhotn/1sz2/gPLHw+m0
aGBE80LKH4/xCz0vnpy5iPLHY+crl3FJ2xdT/nikNvFFpYrGlD8etfr+3tkn
TSh/PO4W1znlh5lS/njsfs1nfhSB8sfD5+q+f3ddNaP88TDOrrBbIF1K+eNh
5VSuUKdhTvnj8WtSxddKkQXlj8fPUF/9l3+toPx82JmXezqeXUn5+ci/ygQM
3WpF+fkQLnTr8dplTfn5kNhUHY18Zkv5+fjSM+Lgf4Ucys+H4JS6xr5hdpSf
j/t9i09NyF5L+fnQ3FRtn2flQPn52OZ+vdZT15Hy87Hi6tp7IRbOlJ+PwC/K
/8za5Ur5+Zj7+hzvxTh3ys/HKLl7Zu+DuZSfj7RZovmjWz0oPx+25v9MKeV7
UX4+Gi/tCxw5yYfyCyCdcaZd554v5RcgK3e/uqevP+UXING1/kXLuwDKL8CY
lXEDovWCKL8Au5hPisnawZRfgPZxZ05PPUXmCiCS99wsOUfmCXBgbnKKmQGZ
ESC5eddn9cX0+zIFGDJ085FDjdRPLMBV8zEquvPo85QIcEih47q2H31eiQBW
WrfT12V5U34BVItTSmPvelJ+ATLPyhkPUKbj0SVA489s3275DZQ/AZZhw1W+
ablR/gRsWlycnr+Njj87AdoOywrkFtL50U9Az4YzQ0w16PwhAZ+qJkVeNqfz
y0lA6eg9Zrn711D+BDTvvTBuFZuuD14C1i1+W6+834byJ+Df58EXDExWUf4E
RI57fVMwyZLyJ0ChZ9l7r4V0vZYkYLZ8QdpaIV3PkgScyun1E31ZRvkT4Hah
4UxzFt0P0gTI2XVvGbma7peuBASsTBxoMJPuJ1YiWreWp1SaLqH8iQjNubt4
hwvdn+xE+O076ZV2iO5f/UQIy3vLjIfS/Y1EbFkYO6A0W3b/JyLE/UFZphnN
D9xEbJ9j7nx6MM0fvETc8Tn32O8dzTdMIizZJ1qcn9N8lJkI94pp1wLf0nwl
TsSUp7s6r/2m+a0kEWk3nv88okOWJCJ6XYDt6oc0P1YnwtF2TmqrO82f0kTU
NdlevDBkFuVPxLSaz6f3VtH8y0pC3FW/6mOHaX5WSkL6AfG9Jek0f7OTUM3r
23RgC83v+kmYVnDRfU/2VMqfBPObqv88OUvPC04SCjxneuq9oOcJNwlpExOu
HFUm85KgmDVFx8KJnkdMEjyv3+6bV0DPq8wk6E3rMTzFIouTsHXn4GvfvGTP
vySUOd+Iiaum56MkCYzNjN6dFuTqJGjWzWgTbqbnqzQJly/tyh2WQc/frj/9
7psXNR+j5zMrGd/fO/YaPqDnuVIyxqQuquT8ouc9Oxkdn4pzGWPZ8z8Z8csG
HZ2dQusDJCOnrOdVVROtJzjJMPS6pPNiEZmbDIFLgLDrMK0/eMnQUf6wJE+F
zCRjm961rdlZtH7JTMZN1bYW0ViyOBnnz6Z1auXL1j/JmJR4fvnSxWRJMth2
SsyZBlovVSdjXvDHhZxksjQZ3OhKGwVdclcyfn4ca/34La2/WAz4H1TlC46Q
FRjcPZK37mAgWYkB9m9werqIrMZg8Pu8tLBRZDaDXblK7ue6af03nYHlzGqj
T01kfQZ6rUeQ8pBsyKDMeMeEH/dl68c/P//gyZj2avJKBkcFTbdHNcvWkwym
1GdMdpSjfk4MDnp8efpolmx9ySBDdHpEpS/Zn8HQW8yANyWy9SaDhCrVe/LD
KX8cg5k6GkuYRDLD4FCk3/BFg+h4pjFQK36lW1pAzmSwyWPSBV4AnY8cBipj
XwY1rKbzJ2bwUNDp+dWDzvfxP8czLuVy/XG6PkoY5F18MmjoHLq+LjMw+aly
s2MwXY8SBl8fTl2zW5Ou3woG744bt7SI6HqvZnBhhooSbzHdLw0M3HVUPA4E
0P0mZfDxxd+q483pfm1l4GboW6CuJbvfGQwU2xSoN9J88R+DzJJRd1M5NB+x
hJih7jpZOJTWGwpCHBDWZCQG0HytJISmi730XhDNz2pC9GBVnP3s5f1mC7Hq
9uqzaQdp/p0uxI+PUcNXaNF8qi/Ew4z8baxImg8NhWALHljPqJGtH4UwaW3q
0Gij47NSiNBdxxe+WUPHgyNE6qHLrEoD+rxOQmjsUlBpqKf1IleIpMgvgwxN
6fj7C5EwcV65S1mt6f/MEyJdTk57vK8q/uc4Ic4FuvZImjX7zQhxQsFxS/m+
Kf1OE6I0qb7viY52vzOFEOi2e/1Xq9bvHCH4noalDzYp9FssxBI1q+waoVK/
jwvhvzW0epfquH6XCGEzc4XxxKPj+31ZiK7W/TH7Xsj2B4W4Nk2x4kozvY9V
CGF7eMeA7yqy9zkhmh+Vt7kNpvfHBiHcND6GfHWm91WpEA7ln82uD6f339Y/
5+/dl93Vg//qd5cQ3kulCplO9D7/nxDxnaMa5etpP4ElwvJhrZuXbqL9CQUR
Zt0f5du1QbYfIkLIHuPWQ+PJaiIYNRkeLbel/Re2CA+3922w1aL9oOkibJav
ddtcQPtR+iJ8CDm80nAN7WcZisDJOF9XNEG2fyaCfEh6RIs87betFGGxXLFP
xADaz+OIkPLG5nzAeNofdBKBcb6quNKC9h+5Iii9ncmemkz7nf4iHN125nZe
FR1fngiTJ9TaOE4lx4kghUVwU7ps/1iEzDcFbzpY5DQROpsORr7ro/3qTBFC
W5NNliqRc0QoZcS52Xq0Hy4W4dePs1Zn3Gl//bgI08adLl2US/v3JSJYfe82
+PmOvh+4LEKDieueVtn3CxIRtj6KivpLTN9HVIjQfqxxLH9Y283/AwpThC4=

      "]]}}, {}, {}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameStyle->Medium,
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 601.}, {-229.6172511018995, 10.802280704378365`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.7846812881123705`*^9, 3.784791771993388*^9, 
  3.7848403358203464`*^9, 3.7848429182908297`*^9},
 CellLabel->"Out[52]=",ExpressionUUID->"4a757256-3a64-41c1-9726-d8a277ee43f5"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"res", "[", 
  RowBox[{"[", "438", "]"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"res", "[", 
  RowBox[{"[", "548", "]"}], "]"}]}], "Input",
 CellChangeTimes->{
  3.784842962118108*^9, {3.7848430367058716`*^9, 3.7848430411141386`*^9}, 
   3.784843271595895*^9},
 CellLabel->"In[56]:=",ExpressionUUID->"4fc8b992-c59f-4380-a0a5-681252dc701f"],

Cell[BoxData[
 RowBox[{"-", "229.6172511018995`"}]], "Output",
 CellChangeTimes->{3.7848404762888317`*^9, 3.784843041660985*^9, 
  3.7848432717208595`*^9},
 CellLabel->"Out[56]=",ExpressionUUID->"b9b453c2-5945-4bd5-a29b-85e4974b02b0"],

Cell[BoxData[
 RowBox[{"-", "39.2385466470192`"}]], "Output",
 CellChangeTimes->{3.7848404762888317`*^9, 3.784843041660985*^9, 
  3.7848432717365103`*^9},
 CellLabel->"Out[57]=",ExpressionUUID->"317ffa6f-36d8-4749-a408-432e3a25c3a7"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Clear", "[", 
  RowBox[{"h0", ",", "\[Phi]"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h1", "=", 
   RowBox[{"2", "\[Pi]", "*", "0.001", "*", 
    RowBox[{"(", 
     RowBox[{"438", "-", "300"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h2", "=", 
   RowBox[{"2", "\[Pi]", "*", "0.001", "*", 
    RowBox[{"(", 
     RowBox[{"548", "-", "300"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"H1", "[", "t_", "]"}], ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"h1", " ", 
       RowBox[{"Sin", "[", "\[Theta]", "]"}]}], "+", 
      RowBox[{"\[Epsilon]", " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Omega]", " ", "t"}], "]"}]}]}], ",", 
     RowBox[{"\[Epsilon]", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"\[Omega]", " ", "t"}], "+", "\[Phi]"}], "]"}]}], ",", 
     RowBox[{"h1", " ", 
      RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"H2", "[", "t_", "]"}], ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"h2", " ", 
       RowBox[{"Sin", "[", "\[Theta]", "]"}]}], "+", 
      RowBox[{"\[Epsilon]", " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Omega]", " ", "t"}], "]"}]}]}], ",", 
     RowBox[{"\[Epsilon]", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"\[Omega]", " ", "t"}], "+", "\[Phi]"}], "]"}]}], ",", 
     RowBox[{"h2", " ", 
      RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LLG1", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mA", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"H1", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mA", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mB", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mAz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mAx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mA", "'"}], "[", "t", "]"}]}]}]}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mB", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"H1", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mB", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mA", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mBz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mBx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mB", "'"}], "[", "t", "]"}]}]}]}]}], "]"}]}], "}"}], "//",
       "Flatten"}], ")"}], "~", "Join", "~", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"mAx", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mAy", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mAz", "[", "0", "]"}], "\[Equal]", 
       SqrtBox[
        RowBox[{"1", "-", 
         RowBox[{"\[Delta]m", "^", "2"}]}]]}], ",", 
      RowBox[{
       RowBox[{"mBx", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mBy", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mBz", "[", "0", "]"}], "\[Equal]", 
       RowBox[{"-", 
        SqrtBox[
         RowBox[{"1", "-", 
          RowBox[{"\[Delta]m", "^", "2"}]}]]}]}]}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LLG2", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mA", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"H2", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mA", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mB", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mAz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mAx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mA", "'"}], "[", "t", "]"}]}]}]}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mB", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"H2", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mB", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mA", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mBz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mBx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mB", "'"}], "[", "t", "]"}]}]}]}]}], "]"}]}], "}"}], "//",
       "Flatten"}], ")"}], "~", "Join", "~", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"mAx", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mAy", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mAz", "[", "0", "]"}], "\[Equal]", 
       SqrtBox[
        RowBox[{"1", "-", 
         RowBox[{"\[Delta]m", "^", "2"}]}]]}], ",", 
      RowBox[{
       RowBox[{"mBx", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mBy", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mBz", "[", "0", "]"}], "\[Equal]", 
       RowBox[{"-", 
        SqrtBox[
         RowBox[{"1", "-", 
          RowBox[{"\[Delta]m", "^", "2"}]}]]}]}]}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"HFM", "=", 
   RowBox[{"ParallelTable", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"\[Phi]", "=", 
       RowBox[{"2", "\[Pi]", "*", "0.01", "*", "j"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"sol", "=", 
       RowBox[{"NDSolveValue", "[", 
        RowBox[{"LLG1", ",", "variables", ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0", ",", "tmax"}], "}"}], ",", 
         RowBox[{"MaxStepSize", "\[Rule]", "0.001"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Is", "=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "2", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "1", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "1", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "2", "]"}], "]"}], ",", "t"}], "]"}]}]}], ")"}], 
        "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "5", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "4", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "4", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "5", "]"}], "]"}], ",", "t"}], "]"}]}]}], 
         ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"NIntegrate", "[", 
       RowBox[{"Is", ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", 
          RowBox[{"tmax", "/", "2"}], ",", "tmax"}], "}"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "0", ",", "100"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"SF", "=", 
   RowBox[{"ParallelTable", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"\[Phi]", "=", 
       RowBox[{"2", "\[Pi]", "*", "0.01", "*", "j"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"sol", "=", 
       RowBox[{"NDSolveValue", "[", 
        RowBox[{"LLG2", ",", "variables", ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0", ",", "tmax"}], "}"}], ",", 
         RowBox[{"MaxStepSize", "\[Rule]", "0.001"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Is", "=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "2", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "1", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "1", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "2", "]"}], "]"}], ",", "t"}], "]"}]}]}], ")"}], 
        "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "5", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "4", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "4", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "5", "]"}], "]"}], ",", "t"}], "]"}]}]}], 
         ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"NIntegrate", "[", 
       RowBox[{"Is", ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", 
          RowBox[{"tmax", "/", "2"}], ",", "tmax"}], "}"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "0", ",", "100"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"ListPlot", "[", 
  RowBox[{"HFM", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Black", ",", "Thick"}], "}"}]}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"ListPlot", "[", 
  RowBox[{"SF", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Magenta", ",", "Dashed"}], "}"}]}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.784667182320609*^9, 3.784667253820386*^9}, {
   3.784676643439802*^9, 3.784676753292672*^9}, {3.784676839261799*^9, 
   3.784676883296608*^9}, {3.784681439703745*^9, 3.7846814472283354`*^9}, 
   3.784682509117309*^9, {3.7846847013259273`*^9, 3.7846847150930715`*^9}, 
   3.784685734479961*^9, {3.784686626161195*^9, 3.7846866387803526`*^9}, {
   3.7847919049766836`*^9, 3.7847919240420675`*^9}, {3.784791977685645*^9, 
   3.7847920005803595`*^9}, {3.784840516338649*^9, 3.784840521267927*^9}, {
   3.784843295387545*^9, 3.7848432994211445`*^9}},
 CellLabel->"In[58]:=",ExpressionUUID->"a354ad9c-c7e5-4ac0-9224-66c264ac14ac"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {GrayLevel[0], PointSize[0.009166666666666668], Thickness[Large], 
     LineBox[CompressedData["
1:eJw11QtMk1cUB/AiiugiYEBAcRs4QGGAgrwEgX+Bgi2vtrRQUHmU54aGDGS8
NHwioIFQwcmEDTTLBjMsQxiDmc1pF1QC4nQ8NgYqZPLckAFR2HhlLj23SdP8
+s/9+vXec85noUyXJm3g8XhJr9//f2pecz51TgUlfJcTavoCWrzSzFoHZl3c
yBp7brqL2QBrAmgPL6WRjZBWrBCsdjKbwq7+dNy1UubdmJk7pm0KZnOM2w/E
XJx6n7wHM9EvJbFFzJY4c32uYsyQ2RquOYl556+8R96HMJ8NKz/oM9tC+3jb
Dq0zqWQ7CPaF1I2OpJAdcPec8y+z7swHEDqqWuu6kEx2xDs5JnM7HyeRnWCa
WPHpRX3mg+icHrEbC0wkO6O2yXF0MDeB7IJXSvsY1wYl2RVTJatfPXgYT3ZD
4OTGnpPzcWR3JGzc3724nfkQeLaetqVxMWQPaLUVZWd9cozsCf0S6+ycvmjy
YdQN163Hbosie0HZKzixLIwke+MpX2hkXCYn++Ds6RLPgt5wMiC99mxlzUKq
MQdMfDcbL8sXU85HiVmPXp46hHI+DA73ND36QkS5L3y6d9QvNARS7ovwMZ3v
jbr9KfdDj6mTwmW7L+V+WFOmt9Sse1Puj96gQwsHHTwo90dHeeW3ae0ulAvw
0ZeV8ya1jpQLcNZscfpo57uUByA+Q9SQvceK8gB4Jfepyz57i/JA/F5b3Rfc
YkR5INzmL8td13UpP4KizxP+nAj7544mPwKDLbmNNem/acwT4orMOmDFbsBH
kwuhskwxyXX6Q2OeCHUva0qK9/5KuQi3JaM/lucP0PogZJTnx22uWqLrB+HW
QGuFRyH7/WDcap2aECay+wvGzJZJp34Bu/8QRHPCAx+8sKQ8BEOqmvL4Mvb/
Q+F+zv3ynTy2P6FQKB89DLzE9i8MBjerXVQ6bH/DXvfzvarN/Wz/xYi1tYop
W+ZrDDHUWoaNZt/Q+XFiGMtcMx6U0/mqxTDr1YWqlJ2/BK2DuZNu9VQfkODv
OftQnKT64STgGQ94FBlSfakl2Gbu9mZTB6s/KdqTHO+v5FN9QorZzKP+9R5U
v5wUIr3Eu8/WFLReisY3lu2GOlj9h2PwrxtuXaXUHwiHiU2IeERM/cOFI7O5
5qrNVuovdTiMZ3tamqdZv8ngU/Ch3Po+9SdkaDO8ntl8lfqXk0HBEyzkZlB/
q2VIXSm93QrW/3KorFY9q3RpPkAOrawlneNdZE4Ou77nFZLCZFovh3Ol8+Of
nNj8iUBV6NrNTUNkRKC9h79cmE3zi4vAqX9f7K7eSvNNHYG944tPOi+xeReJ
4hB9TqVP8xGRWH0S2ZTOkblIeI6c+vrtcbI6EoUfb7LZ5cXmrwKpyXqzlhfI
UOBnkffTV/fInAJWFsPmskWyWgGlX3+o4U42/6Oguz/lvCF7PiAKKqXm+fEf
W+nmwg==
      "]]}}, {}, {}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 101.}, {-225.6642487037647, 0.5686612565124436}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.784667546363476*^9, 3.7846820722165203`*^9, 
  3.784683400794939*^9, 3.784685578515653*^9, 3.7846875409476137`*^9, 
  3.7847926455146275`*^9, 3.7848411623383455`*^9, 3.7848439346786585`*^9},
 CellLabel->"Out[67]=",ExpressionUUID->"4fecd13f-0025-4b6a-8e12-786a3446a910"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {RGBColor[1, 0, 1], PointSize[0.009166666666666668], AbsoluteThickness[
     1.6], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJw11QlM01ccB/Cim0XuASuok1ujCEoBAZHjy1GglKMthRaDsing2GY4ZIsm
HP+EbYEGLAwHEx0GFtCgeI2NaRQ6kMzJJRWPbWKRsakgyMChYAg7+L0mTfP5
f9/v/3/v9b33d9ybJU9fwePx0v/9/ve7/JkOfj2ZZxrpHKSlC7hi577Y8zaz
ISobtthfNGa2gJmttHiRz2yNhRF+VpMZsy1WVVS9+Ysj8zu4Webg6R/A7ICH
x9TXzXczO6G2RD9t8BmzC6b5LpUtbcwbUZk2MN0xxbwJNw4siUVewWRXzAp1
1Vwxsxuqm2Q+CXrmrVAOb/g4ywlkD/xcrm5/rmAWYn2QZdlHA8yemF1d0Lvk
EUL2ws7fra/lnGD2xnFB9rjzmlDydhxp39GVco7ZB4HRH8QspoaRffH35Ix7
/QyzH74433NvXhpO3oFbdxwEL9uY/dG//t17yk0i8k7kFV4dEjQwB0CuNZmK
dY0gB+Kr+oCQXa3MQdCnvWH6IC6SHAzhQJFz4xwzMGSU3fnWvqhlc8C3615s
zJwl80JguKck6vJmas+F4NWzg+uMKln/QuEqaR00OEzzxYUiQvadnafCl/Iw
TBRyZ75OFFIehnN9DTXFY66Uh2PbE73EzXcD5eHo/qZAY9HnQLkIAr2VeJBb
S7kI6uzZWxcuWVMega2lIe4p980pj0BVa7df5EpDyiNxYPhyk+8ZA8oj8dRX
oGmJm+9YzqMwbhbLc//+r2VzUdgzpq0fu/mEcjEyyv+ofvbpKOVi3BXt55UW
j1AejZIf+rxGG1keDU/BpIfB6qeUS2CRfcpuTM/uL0G1Q1rNr/ns+THwz5Xx
+VrWvxiU/1SjEpuw/sfi1KHb+5t0bHyxcDbP0Cg92PjjcPgoX6eXsvmJQ6FB
zuDpk2z+4iF/eKf09mMXyuPxaKHVJErJ5l+KtUVmMWl9HsuGFEbH/Iak3j7U
XgrdmsfcsA/tH60USXqXFUXh7P+Xwa7uQbNuC60PyNAvEVqa5oipXobmo9Oa
qq5oqpchdKrTqfaShOrl+O3DH/W9tmTIYbP7ysvR09Sek6MDm2Ue3XQ/rRz9
qvfCndzJvARs12e2pdbRekUCrt7dO5hly9ZzAmz8WrRetdQ/bQLcFSUz1vZs
/SvgNz7kOmEVQfUKZKpKYjOWaHycAgkN4597WZK1CqQL/xw7G8/2byJWFafm
q8/SfkciHuUNKt0cyVwiCvZ5dzaepPNCm4iKXQJhhT07P5KQm29VN3mN9g+S
cN+4tKfhfbafkjBn21DWbkfWJiHw+sXzXR3sPFPCWneitDeZDCXiutQ2Yj6Z
U2I+Q6PuOk7npVaJzgvBr9qb2PmpwifzkWZTIjJUWGkxdOiFhsypkDKLIzcq
Wb0Kwc3bnpvksvpkfGk8t7jwOpDqkzExcvD/98c/P/3ivg==
      "]]}}, {}, {}, {}, {}},
  
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 101.}, {-55.39532885091877, 0}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.784667546363476*^9, 3.7846820722165203`*^9, 
  3.784683400794939*^9, 3.784685578515653*^9, 3.7846875409476137`*^9, 
  3.7847926455146275`*^9, 3.7848411623383455`*^9, 3.7848439347411556`*^9},
 CellLabel->"Out[68]=",ExpressionUUID->"8f5c1ab7-6e3a-4a25-a6e1-d5e41cd856ed"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"HFM", ",", "SF"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"Black", ",", "Thick"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"Magenta", ",", "Dashed"}], "}"}]}], "}"}]}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", 
   RowBox[{"AspectRatio", "\[Rule]", "1"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.784793068694593*^9, 3.7847930883849087`*^9}},
 CellLabel->"In[69]:=",ExpressionUUID->"65412122-1227-4a6a-9be1-4a806134aec6"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {GrayLevel[0], PointSize[0.009166666666666668], Thickness[Large], 
     LineBox[CompressedData["
1:eJw11QtMk1cUB/AiiugiYEBAcRs4QGGAgrwEgX+Bgi2vtrRQUHmU54aGDGS8
NHwioIFQwcmEDTTLBjMsQxiDmc1pF1QC4nQ8NgYqZPLckAFR2HhlLj23SdP8
+s/9+vXec85noUyXJm3g8XhJr9//f2pecz51TgUlfJcTavoCWrzSzFoHZl3c
yBp7brqL2QBrAmgPL6WRjZBWrBCsdjKbwq7+dNy1UubdmJk7pm0KZnOM2w/E
XJx6n7wHM9EvJbFFzJY4c32uYsyQ2RquOYl556+8R96HMJ8NKz/oM9tC+3jb
Dq0zqWQ7CPaF1I2OpJAdcPec8y+z7swHEDqqWuu6kEx2xDs5JnM7HyeRnWCa
WPHpRX3mg+icHrEbC0wkO6O2yXF0MDeB7IJXSvsY1wYl2RVTJatfPXgYT3ZD
4OTGnpPzcWR3JGzc3724nfkQeLaetqVxMWQPaLUVZWd9cozsCf0S6+ycvmjy
YdQN163Hbosie0HZKzixLIwke+MpX2hkXCYn++Ds6RLPgt5wMiC99mxlzUKq
MQdMfDcbL8sXU85HiVmPXp46hHI+DA73ND36QkS5L3y6d9QvNARS7ovwMZ3v
jbr9KfdDj6mTwmW7L+V+WFOmt9Sse1Puj96gQwsHHTwo90dHeeW3ae0ulAvw
0ZeV8ya1jpQLcNZscfpo57uUByA+Q9SQvceK8gB4Jfepyz57i/JA/F5b3Rfc
YkR5INzmL8td13UpP4KizxP+nAj7544mPwKDLbmNNem/acwT4orMOmDFbsBH
kwuhskwxyXX6Q2OeCHUva0qK9/5KuQi3JaM/lucP0PogZJTnx22uWqLrB+HW
QGuFRyH7/WDcap2aECay+wvGzJZJp34Bu/8QRHPCAx+8sKQ8BEOqmvL4Mvb/
Q+F+zv3ynTy2P6FQKB89DLzE9i8MBjerXVQ6bH/DXvfzvarN/Wz/xYi1tYop
W+ZrDDHUWoaNZt/Q+XFiGMtcMx6U0/mqxTDr1YWqlJ2/BK2DuZNu9VQfkODv
OftQnKT64STgGQ94FBlSfakl2Gbu9mZTB6s/KdqTHO+v5FN9QorZzKP+9R5U
v5wUIr3Eu8/WFLReisY3lu2GOlj9h2PwrxtuXaXUHwiHiU2IeERM/cOFI7O5
5qrNVuovdTiMZ3tamqdZv8ngU/Ch3Po+9SdkaDO8ntl8lfqXk0HBEyzkZlB/
q2VIXSm93QrW/3KorFY9q3RpPkAOrawlneNdZE4Ou77nFZLCZFovh3Ol8+Of
nNj8iUBV6NrNTUNkRKC9h79cmE3zi4vAqX9f7K7eSvNNHYG944tPOi+xeReJ
4hB9TqVP8xGRWH0S2ZTOkblIeI6c+vrtcbI6EoUfb7LZ5cXmrwKpyXqzlhfI
UOBnkffTV/fInAJWFsPmskWyWgGlX3+o4U42/6Oguz/lvCF7PiAKKqXm+fEf
W+nmwg==
      "]]}, 
    {RGBColor[1, 0, 1], PointSize[0.009166666666666668], AbsoluteThickness[
     1.6], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJw11QlM01ccB/Cim0XuASuok1ujCEoBAZHjy1GglKMthRaDsing2GY4ZIsm
HP+EbYEGLAwHEx0GFtCgeI2NaRQ6kMzJJRWPbWKRsakgyMChYAg7+L0mTfP5
f9/v/3/v9b33d9ybJU9fwePx0v/9/ve7/JkOfj2ZZxrpHKSlC7hi577Y8zaz
ISobtthfNGa2gJmttHiRz2yNhRF+VpMZsy1WVVS9+Ysj8zu4Webg6R/A7ICH
x9TXzXczO6G2RD9t8BmzC6b5LpUtbcwbUZk2MN0xxbwJNw4siUVewWRXzAp1
1Vwxsxuqm2Q+CXrmrVAOb/g4ywlkD/xcrm5/rmAWYn2QZdlHA8yemF1d0Lvk
EUL2ws7fra/lnGD2xnFB9rjzmlDydhxp39GVco7ZB4HRH8QspoaRffH35Ix7
/QyzH74433NvXhpO3oFbdxwEL9uY/dG//t17yk0i8k7kFV4dEjQwB0CuNZmK
dY0gB+Kr+oCQXa3MQdCnvWH6IC6SHAzhQJFz4xwzMGSU3fnWvqhlc8C3615s
zJwl80JguKck6vJmas+F4NWzg+uMKln/QuEqaR00OEzzxYUiQvadnafCl/Iw
TBRyZ75OFFIehnN9DTXFY66Uh2PbE73EzXcD5eHo/qZAY9HnQLkIAr2VeJBb
S7kI6uzZWxcuWVMega2lIe4p980pj0BVa7df5EpDyiNxYPhyk+8ZA8oj8dRX
oGmJm+9YzqMwbhbLc//+r2VzUdgzpq0fu/mEcjEyyv+ofvbpKOVi3BXt55UW
j1AejZIf+rxGG1keDU/BpIfB6qeUS2CRfcpuTM/uL0G1Q1rNr/ns+THwz5Xx
+VrWvxiU/1SjEpuw/sfi1KHb+5t0bHyxcDbP0Cg92PjjcPgoX6eXsvmJQ6FB
zuDpk2z+4iF/eKf09mMXyuPxaKHVJErJ5l+KtUVmMWl9HsuGFEbH/Iak3j7U
XgrdmsfcsA/tH60USXqXFUXh7P+Xwa7uQbNuC60PyNAvEVqa5oipXobmo9Oa
qq5oqpchdKrTqfaShOrl+O3DH/W9tmTIYbP7ysvR09Sek6MDm2Ue3XQ/rRz9
qvfCndzJvARs12e2pdbRekUCrt7dO5hly9ZzAmz8WrRetdQ/bQLcFSUz1vZs
/SvgNz7kOmEVQfUKZKpKYjOWaHycAgkN4597WZK1CqQL/xw7G8/2byJWFafm
q8/SfkciHuUNKt0cyVwiCvZ5dzaepPNCm4iKXQJhhT07P5KQm29VN3mN9g+S
cN+4tKfhfbafkjBn21DWbkfWJiHw+sXzXR3sPFPCWneitDeZDCXiutQ2Yj6Z
U2I+Q6PuOk7npVaJzgvBr9qb2PmpwifzkWZTIjJUWGkxdOiFhsypkDKLIzcq
Wb0Kwc3bnpvksvpkfGk8t7jwOpDqkzExcvD/98c/P/3ivg==
      "]]}}, {}, {}, {}, {}},
  
  AspectRatio->1,
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 101.}, {-225.6642487037647, 0.5686612565124436}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.7847930819317613`*^9, 3.784793088806766*^9}, 
   3.7848440158733683`*^9},
 CellLabel->"Out[69]=",ExpressionUUID->"289aa11e-46ee-4803-af4e-72e0dcf34d1d"]
}, Open  ]]
},
WindowSize->{1295, 1037},
WindowMargins->{{110, Automatic}, {Automatic, 0}},
FrontEndVersion->"11.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 5, \
2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 13741, 331, 437, "Input",ExpressionUUID->"5981de4b-bc46-466a-bea8-49eb3958bb2a"],
Cell[CellGroupData[{
Cell[14324, 355, 475, 9, 30, "Input",ExpressionUUID->"a729b156-30d0-405c-ba5f-88b108eccbfe"],
Cell[14802, 366, 9591, 173, 236, "Output",ExpressionUUID->"4a757256-3a64-41c1-9726-d8a277ee43f5"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24430, 544, 365, 8, 52, "Input",ExpressionUUID->"4fc8b992-c59f-4380-a0a5-681252dc701f"],
Cell[24798, 554, 234, 4, 34, "Output",ExpressionUUID->"b9b453c2-5945-4bd5-a29b-85e4974b02b0"],
Cell[25035, 560, 233, 4, 34, "Output",ExpressionUUID->"317ffa6f-36d8-4749-a408-432e3a25c3a7"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25305, 569, 13533, 373, 465, "Input",ExpressionUUID->"a354ad9c-c7e5-4ac0-9224-66c264ac14ac"],
Cell[38841, 944, 2887, 62, 240, "Output",ExpressionUUID->"4fecd13f-0025-4b6a-8e12-786a3446a910"],
Cell[41731, 1008, 2888, 62, 246, "Output",ExpressionUUID->"8f5c1ab7-6e3a-4a25-a6e1-d5e41cd856ed"]
}, Open  ]],
Cell[CellGroupData[{
Cell[44656, 1075, 696, 17, 30, "Input",ExpressionUUID->"65412122-1227-4a6a-9be1-4a806134aec6"],
Cell[45355, 1094, 4335, 89, 367, "Output",ExpressionUUID->"289aa11e-46ee-4803-af4e-72e0dcf34d1d"]
}, Open  ]]
}
]
*)

